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A relation between nuclear forces derived using a phenomenological approach and nuclear 
• efltective field theory (NEFT) is proposed from a renormahzation group point of view. A 

' phenomenological nuclear force (Vph) and an NEFT-based TV-potential (Veft) are satis- 

factorily related to each other through the Wilsonian renormahzation group (WRG) method. 
. It is clearly shown that use of the simple contact interactions in NEFT is adequate to sim- 

^ ' ulate small scale phenomena, and that an NEFT-based A'^A'^-potential (Veft) is free from 

' dependence on the model used to describe small scale phenomena. We study the charac- 

, teristics of Veft from a WRG point of view, emphasizing points that have not previously 

been fully recognized. We also use the Viow k method and a unitary transformation method 
to relate Vph to Veft- It is found that they are not appropriate for this purpose. 

^ ! 

Q ■ §1. Introduction 



The nature of the nuclear force is one of the oldest subjects in nuclear physics. 
In one widely used approach to the nuclear force, one employs the one-pion-exchange 
^ I potential (OPEP) as the well-known long range mechanism and uses a phenomeno- 

^ ■ logical model to describe the short-distance mechanism, which is not well known. We 

refer to this type of nuclear force as a phenomenological nuclear force. At present, 
^ I there are some high-precision phenomenological A^A^-potentials, such as the CD- 

Bonn potential^^ and the Nijmegen potential.^-* 

A different approach to the nuclear force, based on nuclear effective field theory 
(NEFT), was proposed in Weinberg's seminal work^^ and has been studied exten- 
sively.**-* The NEFT approach to the nuclear force is described in detail in several 
references. ^-^ NEFT has attracted much interest because Veft***^ has the following 
desirable formal features that are not possessed by Vph- The first feature, which rep- 
resents the most important aspect of this approach, is that it uses a chiral effective 
Lagrangian and therefore describes a nuclear force in a manner consistent with the 
symmetry properties of QCD, in particular the spontaneously broken chiral symme- 
try. The Lagrangian is composed of effective degrees of freedom (d.o.f.) of hadrons 



*' E-mail: nakkan@rcnp.osaka-u.ac.jp 
**^ We consider NEFT in the form proposed by Weinberg^^ in this work. In another NEFT,''-' 
the nuclear potential is not explicitly derived. We do not discuss this type of NEFT in this work. 

***^ Hereafter, we denote an NEFT-based potential by Veft- Similarly, we denote a phenomeno- 
logical potential by Vph- 
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and is the most general, as long as the assumed set of symmetries exists. When we 
consider low-energy A^^A^-scattering, the nucleon and the pion are the effective d.o.f.; 
heavier d.o.f. (heavy mesons, A, etc.) are integrated out. The second feature is that 
the construction of the nuclear force from the Lagrangian is systematic, following 
a power counting rule. With these two features, NEFT is considered to be model 
independent. Regarding phenomenology, several authors have constructed Veft's 
and have shown their usefulness in reproducing low-energy A'^A^-data.^)"^^ 

Despite the desirable characteristics of NEFT, some naive but yet unanswered 
questions may be raised. One question concerns the use of a contact interaction. In 
NEFT, it is believed that the effect of d.o.f. that are integrated out is absorbed by 
the contact interactions between hadrons representing the effective d.o.f. However, 
the structure of the contact interaction seems to be too simple to simulate physics 
of small scales, when compared to those used in phcnomcnological models. Is the 
use of the contact interaction sufficient to simulate small scale physics? The other 
question concerns the model-independence of NEFT. Why is Veft so special among 
the many A'^A/^-potentials? Is it still one of the many phase-equivalent potentials? 
In fact, the claim that NEFT is model independent is based only on a qualitative 
argument and no quantitative analysis regarding this matter has been made. It is 
very desirable to change this situation. 

The key to answering the above questions is to note the size of the model space 
on which an A''A?"-potential acts. In constructing Veft, one uses a cutoff function 
to restrict the momentum states which the nucleon occupies; Veft acts on a model 
space.*) A typical size of the model space for Veft is considerably smaller than that 
for Vph', the model space for Vph is typically rather large (~ a few GeV) compared to 
the energy region of interest (~ a few hundred MeV) . From a renormalization group 
(RG) point of view, the use of a smaller model space implies a rougher description 
of the system; that is, the system is described by a theory which includes less detail 
regarding small scale physics. Thus we can formulate a scenario for the relation 
between Veft and Vph as follows. There are a number of the VphS which reproduce 
the low-energy A^A^-data. They differ from one another with regard to the ways in 
which they model small scale physics. We start with such Vp^s and integrate out 
the nucleon high-momentum states, thereby reducing the size of their model spaces. 
As the model space becomes smaller, the corresponding potential comes to have less 
information about details of the small scale physics. Therefore, it is expected that all 
VphS eventually converge to a single A/^A''-potential with a sufficiently reduced model 
space. In this way, the model dependence of Vph arising from the treatment of small 
scale physics disappears. Furthermore, the small scale part of that single potential is 
expected to be accurately taken into account by simple contact interactions, because 
the details of the small scale physics are no longer important. After all, Veft consti- 
tutes a parameterization of the single ATA/^-potential. The usefulness of NEFT may 
be assessed by examining how well the NEFT-based parametrization of the single 
A'^A'^-potential works. The single model-space potential is free from dependence on 
the details of the modeling of small scale physics, and therefore its parameterization. 



We refer to a model space with an infinite cutoff as the "full space" . 
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VefT) is also model independent in this sense. This is our scenario for the relation 
between Veft and Vph- If this scenario is shown to be valid, the following conse- 
quences may be realized. First, we can answer the questions posed in the previous 
paragraph and obtain a deeper understanding of NEFT. Second, understanding the 
relation between Vph and Veft-, we can evaluate the role of NEFT in nuclear physics 
from a new point of view. These expected consequences provide a good motivation 
for studying this scenario. 

The purpose of this work is to show that the above-described scenario is indeed 
valid and thereby to propose a relation between Veft and Vph. For this purpose, we 
demonstrate the procedure presented above; we start with some VphS and perform a 
model-space reduction. Then we examine whether the obtained potential is free from 
the model dependence of the Vph and if it is accurately simulated by the NEFT-based 
parameterization, Veft- The remaining problem in the demonstration is to deter- 
mine how we reduce the model space. Some model-space reduction schemes have 
been proposed to this time; the Wilsonian renormalization group method, the 
Viowk method, and a unitary transformation method. However, it has not yet 
been established which should be used in NEFT. We will try each of these methods 
in our demonstration and find the proper one. The criterion for finding the proper 
method is consistency with the basic ideas of NEFT. We examine this consistency 
with respect to two points. First, the proper method should be consistent with the 
method for integrating out the d.o.f. in NEFT; an effective Lagrangian is obtained 
by integrating out the heavier d.o.f. using the path integral method. As a second 
check of consistency, the obtained potential with the reduced model space should be 
accurately simulated by the NEFT-based parameterization. In other words, the ob- 
tained potential should exhibit behavior such that the low-energy constants (LECs: 
coupling constants involved in an effective Lagrangian for NEFT) are natural and 
the NEFT-based perturbation is systematic. Natural LECs and a systematic per- 
turbation (where the meanings of "natural" and "systematic" are specified below) 
are an important basis of NEFT for a convergent perturbation scheme. In fact, this 
criterion can be regarded as a test for examining whether NEFT itself is natural 
and systematic. This is a new type of test for examining the fundamental aspects of 
NEFT. We discuss this point. 

We now describe the organization of the following sections. In ^'2.11 we describe 
the possible methods for the model-space reduction scheme. We also discuss the 
properties of the model-space potentials provided by these methods. In ^2.21 we 
discuss fundamental aspects of NEFT, such as its naturalness. In ^2.31 we give 
explicit expressions for Veft to be used in this work. In ^ we present model-space 
potentials obtained from VphS using the three methods. Comparing these, we find 
the proper model-space reduction scheme. In ^ we discuss the characteristics of 
Veft-, emphasizing points that have not yet been fully recognized. Finally, we give 
a conclusion in ^ 
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§2. Formalism 

2.1. Model-space reduction 

In this section, we will describe the possible methods for the model-space reduc- 
tion. First, we attempt to find a reduction scheme by considering the construction 
of an effective Lagrangian. We start with a Lagrangian Ch in which some d.o.f. \^ 
appear explicitly. (The subscript "i/" means "heavy".) The S-matrix element for a 
given process is obtained from the path integral Z: 

Z = [ V^e'^'^^ ^« . (2-1) 



Now we divide the d.o.f. into "heavy" {^h) and "light" {^l) d.o.f. and suppose 
that we are interested in a system with an energy scale for which only "light" d.o.f. 
are important. In this case, the idea of EFT is to obtain an effective Lagrangian 
written in terms of only "light" d.o.f. Therefore, we integrate out the "heavy" d.o.f., 
obtaining 

Z = [ V^Le'^'^^ , (2-2) 



where Cl {L means "light") is the effective Lagrangian. 

The reduction of the model space for the nucleon states also should be carried 
out following the path integral discussed above. This is the Wilsonian renormaliza- 
tion group (WRG) method. ^'^^ In general, we cannot fully carry out the integration, 
because there are an infinite number of and various types of terms in the Lagrangian. 
Even so, as long as we are concerned with low-energy A^A^-scattering and the A^A^- 
interaction is provided by an A^A^-potential, we can manage to perform the integra- 
tion. In this case, we can start with the path integral given in Eq. H2-l() . where ^ is 
the nucleon field that includes the nucleon momentum states restricted by a cutoff. 
Ah- The Lagrangian Ch includes the starting A^A^-potential, Vh- We separate the 
nucleonic d.o.f. into high-momentum {^h) and low-momentum states {^l) by a cut- 
off Al. Then we integrate out the high-momentum states to arrive at Eq. H2-2() . in 
which the effective Lagrangian Cl includes the low- momentum A^A^-potential, Vl. 
In fact, if we work in the center of mass (CM) A^A^-system, the evolution of the A^A^- 
potential driven by the integration is given by the WRG equation (See Appendix 1X1 
for a derivation.) 

. 1,^.".,.', A^-^V"KA. A) . ,2.3, 

where y^"^) is the A^A^-potential for a given channel (partial wave) a, and M denotes 
the nucleon mass. In V^'^\k',k;p,A), A is a cutoff for the one-nucleon momentum, 
p is an on-shell one-nucleon momentum, p = \/AIE with E being the kinetic en- 
ergy of the two nucleons, and k {k') is a one-nucleon momentum before (after) the 
interaction. Note that we are working in the CM system, and therefore the magni- 
tude of the one-nucleon momentum is the same as that of the relative momentum 
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Fig. 1. For an infinitesimal reduction of the cutoff, the WRG equation represents the evolution of 
the A'^A^'-interaction, as graphically shown above. V {V) is the original (renormalized) interac- 
tion. In the loop, the nucleonic momenta are denoted by fci and fc2, which lie in the momentum 
shell that is integrated out. 



of the two-nucleon system. We use the same notation throughout this work. The 
WRG equation given by Eq. (|2-3() is for a single channel a, but the extension to 
the coupled-channel case is straightforward. Equation (|2-3|) is graphically shown in 
Fig. ^ For an infinitesimal reduction of the cutoff, the interaction V evolves into 
the renormalized one V by absorbing the one-loop graph. In the figure, the loop 
diagram includes the intermediate one-nucleon states oi A — 5A < \q\ < A, where q 
is the momentum of one nucleon. One can obtain Vl corresponding to yl = 71^ by 
solving the WRG equation with the initial condition V = Vh for A = Ah- Actu- 
ally, the solution of the WRG equation, Eq. ()2-3() . is identical to Feshbach's effective 
interaction^^) and is on-shell energy dependent. 

The WRG equation can also be obtained in a different way, as has been done 
by Birse et al}^^ They started with the Lippmann-Schwinger equation for A^A^- 
scattering in a model space with a cutoff A: 

T<-.(.',.;p) = + .(2.4, 

zvr^Jo p'^ — k'^ + le 

Then, they differentiated both sides of the equation with respect to A and imposed 
the constraint that the full off-shell T-matrix is invariant with respect to changes in 
A: i.e., dT/dA = 0. This procedure leads to Eq. ()2-3() . In Birse's procedure, one 
does not have to be careful with regard to the choice of the reference frame in which 
one is working; p, k and k' are the relative momenta of the two-nucleon system, and 
yl is a cutoff for the relative momentum. In case of the CM frame, the meanings of 
the variables in Eq. ()2-3() are the same for the two derivations. In deriving Eq. H2-3|) 
using Birse's procedure, it is essential to use the condition that the full off-shell T- 
matrix is invariant. In fact, this condition is a natural choice if one wants to obtain 
an RG equation which is the same as that obtained from the path integral. Let 
us explain this point. Integrating out the momentum shell in the path integral Z 
does not change the Green functions obtained from Z. In non-relativistic quantum 
mechanics, on the other hand, the Green function for the two-nucleon system is given 

by 

' + ' . V ' . (2.5) 



E- H + ie E-Ho + ie E- Ho + ie E- Ho + ie 

1 1^1 
+ ^ . T- 



E- Ho + ie E- Ho + ie E- Ho + ie' 
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where H (Hq) is the full (kinetic term of the) Hamiltonian and T is the T-matrix. 
If one imposes, for consistency with the path integral case, the condition that the 
Green function is invariant with respect to changes in the cutoff, then this means 
that the full off-shell T-matrix is invariant, as is obvious from Eq. ()2-5() . This is why 
full off-shell T-matrix invariance is a natural choice to realize consistency with the 
path integral method. 

To this point, we have discussed the model-space reduction scheme using the 
WRG equation. In deriving the WRG equation, we started from the path integral so 
that we would perform the model-space reduction in a manner consistent with the 
construction of the effective Lagrangian used in NEFT. In addition to this method, 
two other methods have been proposed for the model-space reductions, the Viowk 
method^^-* and the unitary transformation (UT)^^) method. Although these methods 
are not consistent with the path integral, several authors have argued that these 
methods can be applied to NEFT for model-space reduction. For this reason, we 
regard these methods as possibilities for the model-space reduction scheme in NEFT 
and examine their applicability to NEFT; indeed, it is interesting to consider the 
results obtained using these methods. We briefly describe these methods in the 
following. 

We start with the Viowk method developed by Bogner et al}^^ One representa- 
tion of the Viow k method is the RG equation which enforces the condition that the 
half on-shell T-matrix be independent of the cutoff. The RG equation in this case is 
given 

where the argument of the on-shell energy in V^°'^ is suppressed, because does 
not depend on it if the starting potential is independent of the on-shell energy. The 
T-matrix is denoted by T{k' ,k;p), and the meanings of the arguments here are the 
same as before. A low-momentum potential with a given cutoff can be obtained 
by solving this equation. A practical method to solve this equation is discussed in 
Ref. 12). This method gives a low-momentum potential which is non-hermitian. 

Next, we consider the UT method. We follow the procedure discussed in Ref. 13). 
For a given value of the momentum cutoff, we set up a model space and its com- 
plementary space. Then, we introduce the projection operators onto those spaces 
as 

q\<A, (2-7) 

q\>A. (2-8) 

Now, we derive an effective Hamiltonian acting on only the model space by perform- 
ing the unitary transformation 



V 
A 



(27r)3 



n = u^Hu , 



(2-9) 
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with a condition 

riHX = XHt] = . (2-10) 
The low-momentum potential obtained using the UT method is defined by 

V = T]nr] - r]HoT] . (2-11) 

The low-momentum potential obtained in this way preserves the on-shcll T-matrix 
elements, which has been shown in Ref. 13). (For more details about the deriva- 
tion and features of the UT-based low-momentum potential, we refer the reader to 
Ref. 13).) It is noted that the UT-based low-momentum potential is also obtained 
from the non-hermitian low-momentum potential obtained with the Viowk method 
using the hcrmitization method proposed by Suzuki. -^^^ 

2.2. Naturalness, systematicness and integrability 

In NEFT, we rely on assumptions of naturalness, systematicness, and integra- 
bility. These assumptions are relevant to the size of LECs and are necessary for 
a convergent perturbative calculation following a counting rule. We explain these 
assumptions by considering, for simplicity, a pionless effective field theory (EFT(/)), 
where only the nucleon is dynamical. In EFT(^), all interactions are given by A'^A^ 
contact interactions with 2n (n = 0, 1, 2 ■ ■ ■ ) derivatives. In an S-wave scattering, 
for example, 

Vik', k)=Co + C2(fe2 + + . . . ^ (2.12) 

where C2n is the coupling of the contact interaction with 2n derivatives. All d.o.f. 
other than those of the nucleon have been integrated out and their effects are assumed 
to be accurately simulated by the contact interactions. We refer to this assumption as 
integrability. At leading order, one constructs a nuclear force in terms of the contact 
interaction with no derivatives. At the second leading order, the contact interactions 
with zero and two derivatives are included in the nuclear force. At the m-th leading 
order, the contact interactions with 2n (n = 0, 1, 2 • • • , m — 1) derivatives are taken. 
Because LECs of these contact interactions cannot be determined by the assumed 
symmetries alone, the convergence of the perturbation is not guaranteed from the 
outset. For a convergent perturbation, NEFT assumes naturalness, which we express 
by the condition 

%±i^ « 1 . (2-13) 

The last assumption, systematicness, is defined as follows. Suppose we have a set 

of values of LECs for a A^A^'-potential at a given order. Then, system,aticness means 
that when higher-order terms including new LECs are added to the A^A^-potcntial, 
the values of the LECs for the lower-order terms do not change drastically. In 

*^ In NEFT with the pion integrated out, chiral symmetry is not assumed because the pion mass 
is considered to be very large and chiral symmetry is not a good symmetry in that theory. 
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other words, an A'^A'^-potential to be parameterized by NEFT is assumed to exhibit 
behavior that can be parameterized by a one- (and multi-) pion-exchange potential 
plus a convergent power scries of the nucleon thrcc-momcntum squared, when we 
consider NEFT in which the nucleon and pion are dynamical. In fact, as we will see, 
systematicness is useful criterion to find the proper model-space reduction. 

The assumptions explained above, integrability, naturalness, and systematicness, 
arc not independent ideas; heavier d.o.f. are assumed to be integrable to form an 
interaction that satisfies naturalness and systematicness. Also, naturalness and sys- 
tematicness seem to possess the same meaning. If we consider very high order per- 
turbations, this is probably true. However, in practical cases in which we consider 
a few orders of perturbation, they are not always simultaneously satisfied; we find 
in our results a case in which only naturalness is satisfied, while systematicness is 
not. To this time, many authors have not clearly distinguished these assumptions 
because they are not independent. 

Although these assumptions seem to be reasonable, quantitative tests of their 
likelihood are certainly called for. Even though quantitative studies of NEFT have 
been done extensively during the last decade, we wonder whether these studies are 
sufficiently useful to assess the likelihood of these assumptions. What typically has 
been done in previous studies is as follows. In the case of NEFT, Veft is constructed 
and the LECs are determined at a given order with the use of low-energy experi- 
mental data. In this way, it was found that the size of the LECs is natural. The 
convergence of the perturbation was also studied by going to higher-order perturba- 
tions. It was found that additional terms give smaller contributions at higher orders. 
The higher-order corrections yield better predictions and widen the applicable en- 
ergy region. Although all these results support the assumptions, they only amount 
to examinations of a necessary condition. The problem is that there is no way to 
know whether Veft at a given order accurately simulates the A^A^-potential to be 
parameterized. Even if the naturalness is satisfied at a given order, the situation may 
change when one goes to a much higher order, as is seen in our result given below. In 
this case, the systematicness is not realized. This suggests that the above-mentioned 
examination of (the necessary condition of) the assumptions is insufficient. The 
problem arises from the fact that this bottom-up approach never allows one to know 
the TViV-potential to be parameterized by Veft- By contrast, in our top-down ap- 
proach, we know the A'^iV-potential to be parameterized, which is obtained from Vph 
by integrating out the d.o.f. other than those considered explicitly. Even though we 
do not start from an underlying theory like QCD, we regard the obtained potential 
to be parameterized by Veft for the following reasons: it reproduces the low-energy 
A^A^-data; it is free from model dependence, i.e., dependence on the nature of the 
description of small scale physics, as we see below; its behavior correctly models the 
large scale physics, and the effects of the d.o.f. integrated out. Clearly, our top-down 
approach provides a much better examination of the basic assumptions of NEFT. 
We know how well Veft at a given order simulates the AA-potential to be param- 
eterized, and therefore we can study the likelihood of the systematicness as well as 
that of the naturalness. 
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2.3. NN -potential based on nuclear effective field theory 

As explained in the Introduction, we start with Vph and reduce its model space to 
obtain the corresponding model-space potential, to which we refer as Vm- The model 
space is reduced up to an appropriate size for EFT(7r) (EFT(/)). We use the notation 
EFT(7r) in reference to NEFT in which the nucleon and pion explicitly appear, while 
EFT(/) includes only the nucleon explicitly. We use yEFT{-K){^EFT{f)) to represent 
the iV TV-potential corresponding to EFT(7r)(EFT(/)). The obtained potential, Vm, 
is simulated by Veft{-k) (yEFT{f)) with suitably adjusted LECs. In this subsection, 
we present expressions for Veft(-k) and VEFT{f) that are used in this work. 

For Veft{-k), we use a combination of the OPEP and the contact interactions 
defined in a model space, < |fc| (|fc'|) < /l, as 



The first term is the familiar OPEP, while the rest represent contact interactions. 
The momentum transferred is denoted hy q = k' — k. The quantities gA and f.,^ are 
the axial coupling constant and the pion decay constant, respectively. We do not 
include the Coulomb interaction, because we consider only proton-neutron scattering 
in this work. The LECs of the contact interactions are channel dependent. In this 
work, we are concerned with the ^5*0 and '^Si-'^Di channels. The D2- and D4-terms 
are relevant to the ^Si-^Di channel only. The LECs are A dependent and, in the case 
we use the WRG equation [Eq. ()2-3|) ]. on-shell energy dependent. The expression for 
VEFT{f) is obtained by simply omitting the OPEP from Eq. H2-14|) . 

The A^A^-potential [Eq. (|2-14|) ] we use is not fully consistent with Weinberg's 
counting. Because we consider contact interactions with zero, two, and four deriva- 
tives, we should include more irreducible graphs, such as a TPEP (two-pion exchange 
potential) and an OPEP with more than a single derivative ttNN coupling. Regard- 
ing the contact interactions, the zero- and the two-derivative terms are the most 
general, as long as we are concerned with the ^5*0 and ^Si-^Di channels. However, 
there are other types of four-derivative terms that have not been considered here. 
Nevertheless, we use Eq. (|2-14() for VEFT(7r) to simplify our analysis. In spite of this 
simplification, our result should not be changed essentially, because we employ a 
rather small cutoff value, and therefore the details of the TPEP play only a minor 
role; the effect of the incomplete structure of the four-derivative contact interaction 
is expected to be small. Below we find that this is indeed the case. 

In order to put the iVA^-potential into the Lippmann-Schwinger equation [Eq. (|2-4j) ]. 
it is useful to have expressions of the A^A^-potential for each channel. Such expres- 
sions for the A^A^-potential [Eq. H2-14() ] are presented in Appendix ^ Numerical 
values of the LECs involved are also given in Appendix IbI and in Tabled 



{k'\V\k) 



(2-14) 



with 




(2-15) 
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§3. Result 

This section consists of three subsections. In ^3.11 we test the possible meth- 
ods for the model-space reduction scheme following the procedure described in the 
Introduction. After finding an appropriate model-space reduction, we examine the 
evolution of A^A'"-potential due to the reduction of the model space in i i3.21 We find 
that the model-space potential evolves and comes to have behavior that can be accu- 
rately simulated by contact interactions alone; hence, the pion-exchange potential is 
no longer needed. In ^3.31 we present results for the deuteron properties and discuss 
them. 

3.1. Test of model- space reduction 

3.1.1. Wilsonian renormalization group method 

As explained in the preceding section, we start with a realistic phenomenological 
A^A^-potential (Vph) and then reduce its model space to obtain the corresponding 
model-space potential (Vm) to be simulated by V^px^n) containing suitably adjusted 
LECs. We use three Vphs: the CD-Bonn,^) the Nijmegen I^) and the Reid932) 
potentials. For the proton-neutron ^Sq partial wave scattering, we reduce their model 
spaces following the WRG equation [Eq. up to yl = 200 MeV. The result is 

plotted in Fig. [21 We see that the three potentials, which are originally model 
dependent, are all transformed into essentially the same model-space potential. The 
slight difference between the Reid93 and the others may be attributable to a small 
difference in the low-energy phase shift between them. For this reason, we discuss Vm 








-1.5 I ■ 1 ■ 1 ■ ' ■ 1 

50 100 150 200 

It (MeV) 



Fig. 2. The evolution of various phenomenological A''A''-potentials following the WRG equation. 
The diagonal momentum-space matrix elements of the potentials relevant to the np ^So partial 
wave scattering are shown. The upper three dotted curves represent the bare phenomenological 
potentials: the upper, the Reid93, the Nijmegen I and the CD-Bonn potentials. The lower 
three curves are the model-space potentials (Vm) with A — 200 MeV and p — 10 MeV derived 
from the upper three bare potentials. The solid, dashed and dash-dotted curves result from the 
CD-Bonn, the Nijmegen I and the Reid93 potentials, respectively. The solid and the dashed 
curves are nearly coincident. 
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Fig. 3. Simulation of the WRG-based Vm for the np ^ So partial wave in terms of VEFT(-n)- The solid 
curve represents Vm with A = 200 MeV and p = 10 MeV, and it is derived from the CD-Bonn 
potential following the WRG equation. The dashed curve represents V'i;i?T(7r) i which consists of 
the OPEP and the contact interactions with zero and two derivatives. The dash-dotted curve 
represents Vji'T(^), which consists of Vj^yj-.^) along with the four-derivative contact interaction. 
The couplings of the contact interactions are fixed by fitting them to the solid curve using the 
least square method. The solid and the dash-dotted curves are nearly coincident. 



obtained from the CD-Bonn potential as a representative in the following. Although 
use of the WRG equation introduces an on-shell momentum dependence, the effect 
is negligible for Vm with A = 200 MeV and p ^ 50 MeV. In Fig. El we parametrize 
Vm in terms of Veft{-k)^ namely, the OPEP plus the contact interactions. Here, we 
use ^JpT(^) to represent Veft{tt) including contact interactions with zero and two 
derivatives, while consists of ^Jp'^(^) along with the four-derivative contact 

interaction. We find that V^pj,^^-^ simulates Vm rather accurately. Using V^pj,^^-^ 
yields an almost perfectly accurate simulation. The results are consistent with the 
systematicness assumption of NEFT. The numerical values of the LECs are given 
in Table ^ The numbers on the l.h.s. of Eq. H2-13(l are also shown. From these, 
the naturalness is also found to be realized. The phase shifts obtained with Vm and 
^eft{tt) ^'-'^ ^^^^ channel are in good agreement, as seen in Table ^ This result 
supports the validity of our scenario that Veft and Vph are connected through the 
WRG. 

We continue to reduce the model space up to yl = 70 MeV. First, however, we 
explain why we study a case with such a small cutoff value. One practical reason is 
that, as we will see, different model-space reduction methods lead to quite different 
model-space potentials when we use a small cutoff. For this reason, it is useful to 
study such low-momentum potentials in order to find an appropriate model-space 
reduction scheme. There is also the following formal reason. It has been shown that a 
pionless EFT [that is, EFT(if)] is useful in describing low-energy two-nucleon systems 
in spite of its small model space. *^ Even though V^FTif) results from integrating out 



*•* In studying a system including more than two nucleons, we have to use a model space much 
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Table I. Numerical values for the LECs involved in Veft- The values are fitted to the model- 
space potential Vm using the least square method. Vm is obtained from the CD-Bonn potential 
following the WRG equation. Here, we use, for example, 7.56(— 2) to stand for 7.56 x 10~^. 
For the definition of the LECs in the first column, see Appendix|B| The notation used for the 
Veft's in the third row is explained in the text. 
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-37.1 


IC4/IVC2I 




0.116 
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Cf°i' (fm«) 




-0.254 




-0.255 




8.99 




19.0 



the pion from Veft{tt)j no work has exphcitly shown this. In this work, we address 
this issue by carrying out the following examination. We reduce the model space for 
^eft{it) in the same model space) up to a small size and examine whether 

the obtained Vm is well parameterized by VEFT{f)- We also examine the size of the 
model space appropriate for VEFT{f)- (Supplementary discussion of the derivation 
of VEFT{f) from Vfft{-k) is given in ^4.11 ) Our result is plotted in Fig. 0J We 
adopt p = 10 MeV and 30 MeV as the on-shell momentum. We observe a large on- 
shell momentum dependence. The shift of the potential due to the renormalization 
is larger for p = 30 MeV, which is expected from the WRG equation, Eq. ()2-3() . 



larger than those for Vbft(j^) • Fujii et al. showed that the binding energies of three- and four-nucleon 
systems obtained from exact calculations with a bare A'^A'^-potential cannot be reproduced if one 
uses the corresponding low-momentum A'^A'^-potential for A ~ 400 MeV; even A ~ 400 MeV is too 
small in this case.^^' The situation can be improved by using a larger model space or, alternatively, 
by considering many-body forces generated in a model-space reduction. 
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Table II. The phase shifts obtained from various potentials for np-scattering. The entries in the 
second row are the labels for the potentials. The 'bare' potential used here is the CD-Bonn 
potential, and the model-space potential Vm is obtained from the 'bare' one following the WRG 
equation. The Vbft's are obtained by fitting their LECs to Vm- In the second column, 5 is the 
phase shift, and ei is the mixing parameter for the '^S\-^D\ channel. 



bare 


Vm 


A = 200 MeV 

''eFT{tt) ^EFT{tt) 


Vm 
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48.04 


48.04 
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48.06 








164.47 


164.47 


165.22 


164.46 
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0.00 
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0.00 


0.00 
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0.01 


0.01 


0.01 


0.01 


0.01 


0.02 


0.01 


P = 


30 MeV 




64.42 
137.22 
-0.02 


64.42 
137.22 
-0.02 


66.81 
138.96 
-0.02 


64.34 
137.19 

-0.02 


64.42 
137.22 
-0.10 


64.70 
136.77 
-0.07 


64.42 
137.35 
-0.10 
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k (MeV) 

Fig. 4. The evolution of the WRG-based Vm for the np ^Sq partial wave. The solid curve represents 
Vm with A — 200 MeV. The p = 10 and 30 MeV cases give essentially the same solid curve. The 
dashed curve represents Vm with yl = 70 MeV and p = 10 MeV, while the dash-dotted curve 
corresponds to yl = 70 MeV and p = 30 MeV. 



The parametrization of Vm in terms of ^p'y(^) (the contact interactions with zero 
and two derivatives) is quite good, as shown in Fig. El The values of the LECs 
are hsted in Table Q] We find that the Ojff-diagonal components of Vm can also be 
accurately simulated using the contact interactions whose couplings have been fixed 
by fitting to the diagonal components of Vm- The result is displayed in Fig. El The 
results show that the d.o.f. integrated out, namely the pion and the nucleon high- 
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-3 I • 1 • 1 • 1 1 
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Fig. 5. Simulation of the WRG-based Vm for the np ^So partial wave in terms of the contact 
interaction. The solid curves represent Vm with A = 70 MeV. The upper solid curve corresponds 
to p = 10 MeV, while the lower one to p = 30 MeV. Under each solid curve, there is the 
corresponding dashed curve. The dashed curves represent the contact interactions with zero 
and two derivatives (V^FT(7f))- Their couplings are fixed following the method described in 
Fig.H 








-3 I ■ 1 ■ 1 ■ 1 1 

20 40 60 

Ic (MeV) 

Fig. 6. Simulation of the WRG-based Vm for the np ^ So partial wave in terms of V£prp^^y Plotted 
here are the ojff-diagonal momentum-space matrix elements of the potentials, with kfi^ = 10 
MeV. The other features of this figure are the same as those of Fig. |3 The couplings of the 
contact interactions in V^FT(7f) Eire also the same as those used in Fig. |^ 



momentum states, can be accurately simulated only by the contact interactions with 
natural couplings. We thus find that systematicness is also realized. This means 
that the WRG equation provides an appropriate model-space reduction scheme and, 
simultaneously, that the basic NEFT assumptions discussed in ^2.21 are realized. We 
should recall that use of the WRG equation is also consistent with the integration 
of heavier d.o.f. using the path integral, as we have seen in ^2.11 
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3.1.2. Viowk method 

As in the previous subsection, we start with Vph and reduce its model space, in 
this case using the Vioyjk method. Before considering our result, we briefly discuss 
the results obtained with the Viowk method in previous works. In Ref. 12), it is 
shown in detail that the Viowk method transforms various VphS into essentially a 
single Vm- In the following, therefore, we use only the CD-Bonn potential as the 
starting Vph. A simulation of the Vm obtained using the Viowk method (yl ~ 400 
MeV) using contact interactions is reported in Ref. 18). In that paper, it is shown 
that the shifts of the potential due to the model-space reduction (Vm — V^phV) ^-re 
accurately simulated by the contact interactions. We now discuss our result. As we 
find in Fig. [71 the Va/'s obtained from the WRG and the Viowk method are almost the 
same for A = 400 MeV. Even though the Viow k method has no relation with the path 
integral method for integrating out d.o.f., for this value of the cutoff, this method 
is practically effective as a model-space reduction method in NEFT. We reduce the 
model space up to yl = 200 MeV using the Viow k method and show the resulting Vm 
and its parameterization by Vepx(-k) in Fig. |H1 We find that the Veft{-k) simulates 
Vm rather accurately. For momenta around the cutoff, however, Vm has a large 
curvature, which cannot be simulated with the parametrization used here. 

The situation becomes worse, however, when we reduce the model space fur- 
ther. In Fig. ini we show Vm with yl = 70 MeV obtained using the Viowk method. 
Clearly, the simulation of Vm in terms of the contact interactions does not yield 
consistent behavior. This strongly indicates that the Viowk method is not an appro- 
priate model-space reduction in NEFT. It should be recalled that the Viowk method 
is not consistent with the path integral method of integrating out the heavier d.o.f. 

In case of the UT method, the result is almost the same as that obtained with 
the Viowk method. Therefore, we do not present the result for the UT method here. 



0.5 




-1.5 I ■ 1 ■ 1 ■ ^ ■ 1 

100 200 300 400 

k (MeV) 

Fig. 7. Comparison of the Vm's obtained using the WRG and the Viowk method for the np ^Sq 
partial wave at yl = 400 MeV. The solid curve corresponds to the WRG with p = 30 MeV, 
while the dashed curve corresponds to the Viow k method. In both cases, Vph is the CD-Bonn 
potential. 
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Fig. 8. Simulation of Vm obtained using the Viowk method for the np ^So partial wave in terms of 
Veft{tt)- The solid curve represents Vm derived from the CD-Bonn potential using the Viowk 
method with A = 200 MeV. The dotted curve, shown for a comparison, is the same as the solid 
curve in Fig. |21 The other features are the same as in Fig. |3 




20 40 60 

k (MeV) 



Fig. 9. Vm obtained using the Viomk method and an energy independent contact interaction for 
the np ^ So partial wave. The solid curve represents Vm with /I = 70 MeV derived with the Viom k 
method. The dashed curve represents the contact interactions with zero and two derivatives. 
Their couplings are fixed by fitting them to low-energy phase shift data. 

3.2. Evolution of N N -potential 

Through the investigation carried out in the preceding subsection, we find that 
the WRG equation hkely provides an appropriate model-space reduction scheme. In 
this subsection, we use the WRG equation to further examine the evolution of an 
A^A^-potential and thereby confirm the scenario connecting Vph to Veft and study 
the role played by the WRG equation as the method facilitating the connection. 
Recall that we studied the ^Sq channel in the preceding subsection. Here we study 
the ^Si-^Di channel. It is interesting to study this channel, which contains a richer 
variety of phenomena than does ^Sq; in particular, the pion plays an important 
role. This is a good place to examine the evolution of the potential and to find 
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an appropriate size of the model space for which the pion-exchange potential is no 
longer needed. 

First, we examine the evolution of various VphS by reducing their model spaces. 
In the ^Si-^Di channel, all the Vp^s used here are phase-shift equivalent. Regard- 
ing the deuteron D-state probability {Pd), however, there is a model dependence: 
PD(CD-Bonn) = 4.85%, PD(Nij I) = 5.66%^) and PD(Reid93) = 5.70%.^) It is 
interesting to examine whether all of these A^A^-potentials evolve into essentially a 
single model-space potential, as in the ^Sq case, in spite of this model dependence 
of Pz). 

Plots of the evolution of the CD-Bonn and the Nijmegen I potentials for {'^Di\V\'^Si) 
are presented in Fig. UHl where we use A = 200 MeV and p = W MeV. The OPEP 
is also shown there. We do not plot the evolution of the Reid93 potential, in order 
to make the figure clearer. Including the Reid93 potential would not change the 
discussion below. As we see in the figure, the initial model dependence is small for 
low-momentum components because of the dominance of the OPEP tensor force. It 
is noted that there is a strong model dependence in momentum components much 
larger than those shown in the figure. The shift of the potential due to the model- 
space reduction is not large for our choice of A and p, and the degree of the model 
dependence remains small. Therefore we use the Vm obtained from the CD-Bonn 
potential as a representative in the following. We also observe from the figure that 
the shape of the original and the model-space potentials are largely governed by the 
OPEP. Obviously, we need the OPEP to parameterize Vm for A = 200 MeV. The 
simulation of Vm in terms of VEFT(7r) is quite accurate, as seen in Fig. ^2 The values 
of the LECs for this VEFT{Tr) listed in Tabled while the phase shifts are listed 
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Fig. 10. The evolution of np {^Di\V\^Si) for various Vph^s under the WRG equation. The diagonal 
momentum-space matrix elements are shown. The lower two dotted curves represent the bare 
Vph's with the upper (lower) one the CD-Bonn (the Nijmegen I) potential. The upper two curves 
represent Va/'s with A = 200 MeV and p = 10 MeV derived from these two Vph^s. The solid 
and the dashed curves are the Vm's obtained from the CD-Bonn and the Nijmegen I potentials, 
respectively. The OPEP is also shown by the dash-dotted curve. Note that the scale of the 
vertical axis is different from that in Fig. |5|by an order of magnitude. 
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Fig. 11. The simulation of Vm for np Di\V\^ Si) in terms of VeftCit)- The solid curve represents 
Vm with A = 200 MeV and p = 10 MeV. It is derived from the CD-Bonn potential using the 
WRG equation. For the p = 30 MeV case, the corresponding Vm is indistinguishable in this 
graph from the solid curve. The dashed curve represents Vj^FT(7r)i while the dash-dotted curve 
represents • For the notation, see the caption of Fig. |3 

in Table nil 

In order to simulate a Vm without the OPEP, what is an appropriate value of 
A? Because the OPEP can be expanded in terms of the contact interactions as 

jis^ + "^^) ~ I'^'ir — Q'^/iT^t + ■ ■ ■ ) the expansion is convergent for A < 
For a rapid convergence, a much smaller cutoff value is necessary. Therefore, we 
reduce the model space to have the potential with A = 70 MeV ~ m-^l^. The result 
is shown in Fig. for p = 10 and 30 MeV. We find that there is a large on-shell 
energy dependence between the p = 10 and 30 MeV cases, while no difference can be 
discerned in the graph for A = 200 MeV. A simulation of Vm in terms of the D2^^- 
contact term alone [V^prp^^y, see Eq. (llj-3() for the definition of the D^^^-term] is 

also shown in the figure. Inclusion of the Z^^'^^^-term (V^prp^^^) makes the simulation 
sufficiently good that no discrepancy can be seen in this graph. The simulation is, 
even with this relatively large cutoff value, significantly better than what we naively 
expected from the above considerations. We can explain this as follows. As the model 
space is reduced, the potential becomes so strongly renormalized that the bare OPEP 
plays only a minor role; that is, the bare OPEP contribution is hidden by the shift 
of the potential due to the renormalization. The shift of the potential has a shape 
suitable for a contact interaction expansion. This is the reason that the expansion 
is effective for a relatively large cutoff value. This accurate simulation supports 
the systematicness assumption and the values of the LECs in Table Q indicate that 
naturalness is also realized. 

It is interesting to examine the evolution of driven by the model- 

space reduction, because much more significant renormalization of the potential is 
expected on the basis of the following speculation. The tensor force is known to 
excite a low-energy ^Si state to a high-energy ^Di state, and vice versa. Therefore, 
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Fig. 12. Evolution of A''A''-potential under the WRG equation for np Di\V\^ S\) and the simula- 
tion of Vm in terms of VeftHc)- The solid curves represent Vm's with A = IQ MeV; the upper 
solid curve corresponds to the p = 10 MeV case, while the lower one corresponds to p = 30 
MeV. Each of the solid curves is simulated by the Dj^^'-term alone (Vj^i?T(7f))i ^.s shown by the 
accompanying dashed curve. By including the _D4''^'-term {Vi/pT^f))^ the simulation of the solid 
curves cannot be distinguished from Vm in this graph. For comparison, Vm with A = 200 MeV 
is also plotted by the dash-dotted curve. This curve is the same as the solid curve in Fig. 1111 
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Fig. 13. The evolution of np 1^ Si\V^ S{) for various V^^'s following the WRG equation. The upper 
three dotted curves represent the bare V^h's corresponding to (from top to bottom) the Reid93, 
the Nijmegen I and the CD-Bonn potentials. The lower three curves represent Vm's with A = 
200 MeV and p = 10 MeV derived from the upper three bare potentials. The solid, dashed and 
dash-dotted curves represent Vjvf's derived from the CD-Bonn, the Nijmegen I and the Reid93 
potentials, respectively. 

a sequence of transitions like ^5*1 — > ^Si is renormalized into the central force 

of the ^Si — > ^Si transition after integrating out the high-momentum states. Below 
we show that this is indeed the case. We start with three potentials (CD-Bonn, 
Nijmegen I and Reid93) for (^S'i|y|^S'i). In this case, the initial model dependence 
is clear for the low-momentum matrix elements, as seen in Fig. 1131 We see from 
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Fig. 14. Evolution of np {'^ Si\V\^ Si) under the WRG equation. The solid curve represents the 
original CD-Bonn potential while the dashed (dotted) curve is the corresponding Vm with A = 
200 MeV and p = 10 MeV (p = 30 MeV). 

the figure that they evolve into essentiahy the same Vm- From this point, we use the 
Vm from the CD-Bonn potential as a representative. Our expectation of a strong 
renormalization is confirmed by comparing Fig. ^1 with Fig. 1101 Because of this 
strong renormalization, the on-shell energy dependence of Vm is slightly discernible 
even at yl = 200 MeV (Fig. I14j) . We do not show the result but mention that this 
Vm is accurately simulated by VEFT{Tr) at the same level found in Fig. 01 When we 
reduce the model space up to yl = 70 MeV, we observe a significant evolution, as seen 
in Fig. 1151 Even with this evolution, we confirm again that the naturalness and the 
systematicness are maintained, as seen from Fig. llSl and Tabled Although we do not 
present the result, we mention that the contact interactions determined by fitting 
to the diagonal components also accurately reproduce the off-diagonal momentum 
components. We do not present the result for the evolution of {^Di\V\'^Di) either, 
but we find that the trend of the result is the same as that in the Di\V\^ Si) case. 

As we have shown, the large model dependence of the original Vp^s for the 
^Si-^Di channel have disappeared after the model-space reduction under the WRG 
equation, in spite of the differences in Pd among the original Vp^s. This is not 
surprising, however, because these differences in Pd are due solely to modeling of 
small scale phenomena. This point will be discussed in the next subsection. We 
have also shown that Veft{tt) (^eft{-^)) with suitably adjusted LECs accurately 
simulates the low-momentum Vm obtained using the WRG equation and reproduces 
the phase-shifts, as shown in Table ITII 

3.3. Deuteron in nuclear effective field theory 

We now examine how well Vfft{-k) (yEFT{f)) reproduces the wave function and 
the properties of the deuteron. For this purpose, we derive Vm from Vph using 
the WRG equation, where the on-shell energy is set to the deuteron binding energy 
(B.E.). We obtain Veft{-k) {^EFT{i)) by fitting the LECs involved to Vm- We solved 
the model-space deuteron eigenvalue problem in the momentum space with the use 
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Fig. 15. Evolution of Vm under the WRG equation for np Si\V\^ Si) and the simulation of Vm 
with VEFT^f)- The solid curves represent the model-space potentials with A = 7G MeV; the 
upper solid curve corresponds to the p = 10 MeV case, while the lower one corresponds to p = 
30 MeV. Each of the solid curves is simulated using the contact interactions with zero and two 
derivatives (Vj^FT(jf))i ^ shown by the accompanying dashed curve; each of the dashed curves 
is almost completely indistinguishable from the corresponding solid curve in this graph. For 
comparison, Vm with A = 200 MeV is also plotted by the dash-dotted curve. This curve is the 
same as the dashed curve in Fig. 1141 Note that the dashed and the dotted curves in Fig. 1141 are 
indistinguishable with the scale used in this figure. 

Table III. The deuteron binding energy (B.E.) and the _D-state probability {Pd) obtained from 
various potentials. In the second column, appear the VphS from which the Fa/'s are obtained 
using the WRG equation. The numbers in the third column are not from the original paper but 
from our numerical calculation. We use the extended definition of Po [Eq. I3TI |1 here. 
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B.E. 

(MeV) 


CD-Bonn 


2.224 


2.224 
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2.225 
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2.223 


Nij I 

Reid93 


2.226 
2.225 


2.226 
2.224 
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2.224 
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Pd 
(%) 


CD-Bonn 


4.85 


1.44 


1.55 


1.45 


0.09 


0.07 


0.09 


Nij I 

Reid93 


5.68 
5.70 


1.47 
1.45 


1.58 
1.56 


1.47 
1.44 


0.09 
0.08 


0.07 
0.07 


0.09 
0.08 



of Vm {A = 200, 70 MeV), VEFT(n) = 200 MeV) and Veft^) = 70 MeV). 
The deuteron B.E. obtained in each case is hsted in Table HTll We find that yEFT(n) 
and V^prp^^^ reproduce the B.E. quite well. 

How about the other deuteron properties which are often discussed in the liter- 
ature? They are, for example, the asymptotic S'-wave normalization As, the D/S- 
ratio 7] and the -D-state probability Pd- Actually, a calculation of these quantities is 
beyond the capability of a theory like NEFT, in which a relatively small model space 
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is used. In order to determine these quantities, we need information about details of 
the smah scale physics, which have been integrated out in NEFT. Allow us to explain 
this point further. In order to calculate these quantities, the high-momentum com- 
ponents of the deuteron wave function are necessary. The WRG equation provides 
a relation between the wave functions of the full space and the model space; that is, 
the WRG equation does not change the off-shell T-matrix, and therefore the model- 
space wave function is obtained from the corresponding full-space one by just cutting 
off the momentum components higher than a given cutoff. Thus, no information is 
available from NEFT about the high-momentum components of the wave function. 
This is why some deuteron properties cannot be obtained in NEFT. Of course, if 
we use a sufficiently large model space, the quantities can be calculated with a good 
approximation. However, as found from Fig. 1161 such a good approximation requires 
yl ^ 1 GeV, which is too large for NEFT including only the nucleon and the pion 
explicitly. 

Although many deuteron properties are not described in NEFT, we consider the 
deuteron wave function in a model space to see how well Veft works; the result is 
given in Figs. I17H2U1 In fact, the normalization of the deuteron wave function also 
requires information about details of the small scale physics; NEFT alone cannot 
normalize the wave function. Here, we normalized the model-space deuteron wave 
function such that the normalized S'-wave and the wave function obtained with Vph 
give the same amplitude at A; = yl. We confirmed that the wave function obtained 
with Vm is the same as that obtained with the corresponding Vph (the CD-Bonn 
potential in this case), as should be the case. From the figures, we see that V£;FT(7r) 
and VEFT{f) are successful in generating the wave function. 

Lastly, we discuss the deuteron D-state probability, Pd- As we have seen, various 
VphS evolve into a single Vm under the WRG equation. Therefore, if we extend the 



0.1 



I 0.05 

I 





500 1000 

k (MeV) 

Fig. 16. The radial part of the deuteron wave function in the momentum space. The solid curve 
represents the S'-wave, while the dashed curve represents the D-wave. The vertical axis corre- 
sponds to the radial wave function multiplied by the momentum k. 
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Fig. 17. The radial part ol the deuteron 5*- 
wave functions in the momentum space. 
The solid, dashed and dotted curves are 



obtained with Vm, Vi^^j^j and Vj^S-TCvr)! 
respectively. Regarding the normalization 
of the wave functions, see the text. 




20 40 
k (MeV) 

Fig. 19. The radial part of the deuteron S- 
wave functions in the momentum space. 
The solid, dashed and dotted curves are 
obtained with Vm, V|ft(#) and VEFT(f)^ 
respectively. 
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Fig. 18. The radial part of the deuteron D- 
wave functions in the momentum space. 
The other features are the same as those 
in Fig.[T7| 
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Fig. 20. The radial part of the deuteron D- 
wave functions in the momentum space. 
The other features are the same as those 

in Fig. rrai 



definition of Pd to the model-space case using 

Pd = 3.^ 

the model dependence of Pd disappears after the model-space reduction. In Table 
mil we list the extended Pd's for Vivf's originally from various VphS. Actually, 
this result is expected by observing the full-space deuteron wave function for the 
following reasons: (1) only weak model dependence on the choice of Vph is found in 
the low-momentum components of the deuteron wave function; (2) the model-space 
deuteron wave function is obtained from the full-space one by simply cutting off the 
momentum components higher than A. We can conclude that the model dependence 
of Po comes from the details of the modeling of the small scale physics. Coraggio 
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et al. also discussed the point that the model dependence of Pd disappears after 
integrating out the high-momentum states with the use of the Viowk method. ^^-^ 

§4. Discussion 

All the results presented in the previous section support our scenario, discussed 
in the Introduction, for understanding Veft and Vph in a unified manner with the 
help of the WRG. In the following, we take the viewpoint offered by this scenario 
to understand the relation between Veft and Vph. In H4.11 we discuss, from this 
viewpoint, the relation between EFT(7r) and EFT(^), starting from the chiral effec- 
tive Lagrangian for the pion and the nucleon. We also note that the WRG equation 
indicates properties which Veft should possess. We refer to a Veft possessing such 
properties as a proper Veft- However, some of these properties have not yet been 
fully recognized. In ^4.21 we discuss the expected nature of Veft based on consid- 
eration of the WRG equation. From the same viewpoint, we also discuss the nature 
of a previously constructed Veft in ^4.31 We refer to such a Veft as a conventional 
Veft- Finally, in ^4.41 we discuss what we can learn about a description of a low- 
energy two-nucleon system based on a traditional model approach and that based 
on NEFT from the point of view developed in this work. 

4.1. From EFT(tt) to EFT{f) 

In EFT(7r), we start with the path integral 

Z = j VNVN'^VU e'^ '^^ , (4-1) 

where is the nucleon field, U contains the pion field, and C^^ is an effective chi- 
ral Lagrangian containing the nucleon and the pion fields. The usual procedure 
employed in NEFT is to specify a set of irreducible graphs from following a 
counting rule and to regard the set of the graphs as the A^A^-potential or the nu- 
clear operator. One puts these operators into the Schrodinger equation and solves 
it to obtain physical observables. Because one regards a set of irreducible graphs 
as the nuclear operator, the theory is no longer a quantum field theory but, rather, 
non-relativistic quantum mechanics. Here, we describe the above standard NEFT 
procedure in the language of quantum field theory. 

The standard procedure to calculate the A^A^-scattering amplitude in quantum 
field theory is as follows: (1) calculate the nucleon four-point Green function, G4, 
using the path integral of Eq. ()4-l() : (2) apply the LSZ reduction formula to G4 
to obtain the S-matrix of the A^A^-scattering. This procedure is not equivalent to 
that used in NEFT. With the following modification, the two procedures become 
equivalent. In calculating G4 perturbatively with the use of one uses a "counting 
rule"; that is, one retains only the ladders of a set of irreducible graphs. The set 
of irreducible graphs from is, of course, what the counting rule specifies in the 
standard NEFT procedure. We use G4 to represent the nucleon four-point Green 
function to which the counting rule has been applied. By applying the LSZ reduction 
formula to G4, we can obtain the S-matrix for the A^A^-scattering. In this way, we 



Nuclear Forces from Viewpoint of Renormalization Group 



25 



explained the standard NEFT procedure in the language of quantum field theory. 

Now we consider how to obtain EFT(/) from EFT(7r). If we restrict ourselves 
to a two-nucleon system and need to calculate G4 rather than G4, it is equivalent to 
using the path integral 



with the structure oi Cn given in Eqs. HA-2|) - HA-4|) . The matrix element of the 
interaction Lagrangian contained in between free two-nucleon states is, up to an 
overall sign, nothing more than the set of irreducible graphs which the counting rule 
specifies. Somehow, the pion has been integrated out. Once a path integral of the 
form of Eq. (|4-2jl is obtained, the nucleon high-momentum states can be reduced up 
to a size appropriate for EFT(/) following the procedure discussed in Appendix 1X1 
Then one obtains the effective Lagrangian to be directly compared with the pionless 
Lagrangian; the LECs involved in EFT(if) can be fixed by the comparison. In this 
way, one obtains EFT(/) starting from an effective chiral Lagrangian of EFT(7r). 

It is noted that we do not perform a rigorous procedure, which is far from trivial, 
to obtain Eq. H4-2() from Eq. H4T() by integrating out the pion. Nevertheless, as long 
as we are concerned with G4, our procedure does make sense. 

4.2. Characteristics of N N -potential properly based on effective field theory 

First, we discuss the on-shell energy dependence of Veft- As we have seen in 
the result, Veft has the on-shell energy dependence. The dependence is strong, in 
particular, in the case of a small value of the cutoff. This is a natural consequence 
of integrating out the nucleon high-momentum states by use of the WRG equation. 
In spite of this fact, the procedure employed in previous works to construct Veft is 
as follows. The irreducible graphs from a Lagrangian up to a given order are simply 
multiplied by a cutoff function, and this is used as Veft- The high-momentum 
states which are not considered explicitly in this way are assumed to be absorbed 
by on-shell energy independent contact interactions. Possible energy dependence is 
considered to be eliminated by using the equation of motion. However, as discussed 
by Birse et aZ.,^^) the WRG equation indicates that the on-shell energy dependence 
introduced by the WRG equation is not eliminated by the equation of motion. This 
can also be seen clearly from the result {e.g., see Fig. |^ in the previous section. 
Birse et al. treated the on-shell energy dependence coming from the nucleon high- 
momentum states that have been integrated out. However, they were interested in 
the behavior of the potential around the fixed point, and therefore they did not study 
how important the on-shell energy dependence of Veft is- 

We showed that Veft is largely on-shell energy dependent, in particular, in the 
case of a small value of the cutoff. Actually, the large on-shell energy dependence 
is due to the fact that separation of scale is not so strong. In order for contact 
interactions to accurately simulate the pion contribution, we have to take the cutoff 
value to be maximally yl ~ 70 MeV. On the other hand, we hope EFT(^) to be 
effective for p ^ 70 MeV. Therefore, we find p <j A. The WRG equation [Eq. (|2-3|) ] 
immediately indicates that there must be a large on-shell energy dependence in such 
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a case. By contrast, the WRG equation tells us that the on-shell energy dependence 
is negligible in the case p <^ A. 

There is another source of on-shell energy dependence of Veft- The on-shell 
energy dependence due to the recoil correction has been considered in a previous 
work.^) This effect is a higher-order correction and, as discussed in Ref. 7), it may 
be eliminated when one uses a unitary transformation method in constructing Veft- 
This type of on-shell energy dependence is not as important as that from the model- 
space reduction. 

Now we know that it is important to consider the on-shell energy dependence 
of the couplings of the contact interactions in Veft- If we also parameterize the on- 
shell energy dependence, we have to modify the counting rule. A parameterization 
for a pionless on-shell energy dependent A^A^-potential is 

V{k\ k- p, A) = Ci, + {k^ + k'^) +C^2P^ + --- , (4-3) 

as given in Eq. (1) of Ref. 11). A natural choice is to use a modified counting rule 
which regards p2rap^2?ife ^^^j^ p2mpf^2mf, ^gj-j^g ^^^^ _|_ = j^j^ _j_ ^^^^ g^g same order. 

In order for the modified counting rule to be effective, it is necessary that the on-shell 
energy dependence of the LECs can be expanded in a convergent power series in p^. 
In fact, such an expansion is expected to exist. This is because the WRG equation 
[Eq. (|2-3|) ] indicates that the shift of the potential due to a change of the cutoff has 
an on-shell energy dependence which can be expressed as a convergent expansion 
in terms of p'^, if the starting potential has such a dependence. Therefore, if we 
start with a potential without on-shell energy dependence and reduce the model 
space according to the WRG equation, the obtained model-space potential would 
have such a tractable energy dependence. Indeed, we can see from Figs. OTIandl^ 
that the LECs have such on-shell energy dependence. In the figures, we plot the 
on-shell energy dependence of the LECs contained in V^pj,^^^ {A = 200 MeV) for 
np ^5o-scattering; the definition of the LECs is given in Eq. (jBTI) . The solid curves 
are obtained by fitting the LECs to Vm for each p. We simulate the jj^-dependence 
of the LECs using 

C^f^' ip') = + C^^P" + Ct^l^p" + • • • . (4-4) 

In order to obtain the dashed (dotted) curve in the figures, we took the first two 
(three) terms in the r.h.s. of Eq. 1)4 -41) and fitted their parameters to the solid curve 
using a least square method in the region < p < 150 MeV. We restrict the p- 
region in the simulation for the following reasons: (1) the LECs suddenly change as 
functions of p near p ^ A = 200 MeV, and the higher-order terms in Eq. (|4-4|) are 
necessary for a good simulation; (2) inclusion of the high-p region in the simulation 
makes the simulation worse even in the low-p region. 

Next we discuss the model independence of NEFT. As mentioned in the In- 
troduction, the model independence of NEFT is based on the following two facts: 
(1) the most general Lagrangian consistent with the assumed symmetries (which 
consists of effective d.o.f.) is used; (2) a systematic and perturbative calculation is 
performed following a counting rule, and naturalness, systematicness and integrahil- 
ity are necessary assumptions to organize the perturbation. Although the validity 
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of these assumptions had not previously been fully confirmed, we found that they 
are indeed realized in the proper Veft- Therefore, our finding supports the model 
independence of NEFT based on the above-mentioned two facts. However, we have 
the following question: Although NEFT is a model-independent framework in this 
sense, is Veft still one of many phase-shift equivalent potentials defined in the same 
model space? The answer is 'yes'. As we have seen, the model dependence of a 
potential with a sufficiently reduced model space does not arise from the descrip- 
tion of small scale phenomena. Instead, it arises from the choice of the model-space 
reduction scheme, and this is why the answer is 'y^s'- However, we have no choice 
of the reduction scheme if we work within NEFT. In order to maintain consistency 
with the basic ingredients of NEFT, we must choose the WRG method. With this 
renormalization scheme, there is an essentially unique model-space potential to be 
simulated by Veft- 

Lastly, we emphasize one point: We should be careful that Veft is defined in 
a model space with a relatively small cutoff. Of course, some aspects regarding this 
point have been treated in previous works. For example, the cutoff dependence of 
NEFT-based predictions of observables have been often examined. However, there 
has been some confusion resulting from an insufficient understanding of this point. 

One point of confusion is seen in a study of resonance saturation. In Ref. 20), the 
authors started with a phenomenological one-boson-exchange potential and directly 
expanded the one-boson exchange terms, excluding the OPEP, in terms of a series 
of contact interactions and compared the couplings obtained in this way with those 
obtained using the conventional procedure of EFT(7r). However, this comparison 
makes no sense. The potential obtained in this naive expansion cannot reproduce the 
scattering observables generated by the original potential. It should be noticed that 
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Vph and Veft are defined in different model spaces and that tlierefore tliey cannot 
be compared with each other directly. In contrast to the above-stated manner, it is 
necessary to perform the model-space reduction before the expansion. 

Another point of confusion is seen in the determination of the couplings of the 
contact interactions, which is relevant to the ^Si-'^Di partial wave scattering, with 
the use of the deuteron properties. Some authors have used the deuteron binding 
energy (B.E.), the asymptotic S'-wave normalization (As), and the L'/5'-ratio (77) 
to determine the couplings. However, as we discussed in the preceding section. 
As and r] include information concerning the details of small scale physics, and 
therefore they are not well-defined quantities in the model-space NEFT framework. 
We should recognize that the use of As and t] to determine the LECs means that 
we must introduce a model into the NEFT framework. 

4.3. The conventional Veft 

We now discuss the nature of the conventional Veft constructed in the previous 
work from our point of view. The procedure to derive the conventional potential 
is described in the Introduction. We describe their characteristics in the following. 
They are hermitian, energy independent (except for possible higher-order correc- 
tions) and defined in a model space with yl < 1 GeV. The unknown parameters 
they contain are determined such that physical observables are accurately predicted. 
Specifically, this means that the on-shell T-matrix elements are accurately predicted. 
The potential generates an orthogonal set of wave functions. The above-stated char- 
acteristics are the same as those of a model-space potential obtained using the UT 
method. Therefore, the conventional Veft can be regarded as a parametrization 
of a model-space potential obtained using the UT method. Is there another inter- 
pretation? According to the basic idea of EFT, two NEFTs defined in different 
model spaces should be connected by integrating out heavier d.o.f. If one wishes 
to reduce the model space while keeping the conventional Veft equipped with the 
above-mentioned characteristics, the only known method is the UT method. There- 
fore, we adopt the interpretation that the conventional Veft is a parametrization of 
the model-space potential due to the UT method. 

With the above interpretation, we find the following problems in the conven- 
tional Veft- Even if this interpretation is not correct, a considerable number of the 
problems would remain. The proper Veft discussed in the preceding subsection, by 
contrast, is free from these problems. We now consider these problems. In order 
to see the problems clearly, we often consider cases with rather small model spaces. 
First, as discussed in ^2.11 we note that the UT method is unrelated to the path 
integral method of integrating out heavier d.o.f., while an effective Lagrangian for 
NEFT is considered to be obtained using the path integral method. 

Secondly, as discussed in ^ the UT method gives a model-space potential which 
cannot be simulated accurately by the conventional NEFT-based parametrization, 
in particular, in the case of a small cutoff. However, this does not mean that a 
conventional on-shell energy independent Veft cannot predict physical observables. 
For example, suppose that we construct a conventional V^prp^^y for which we have 
two contact interactions with zero and two derivatives, respectively. In this case, 
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we can find values of tlie LECs of the contact interactions which give the scattering 
length and the effective range. As seen in Fig. however, such a potential (the 
dashed curve) cannot simulate the potential obtained with the UT method (the 
solid curve),*) and it is applicable to only a very limited energy range that is far 
below the cutoff. Although we have freedom in choosing the unitary transformation, 
the small applicable energy range indicates that no unitary transformation gives a 
potential that can be accurately simulated by the two contact interactions. This 
is an example of a case in which the naturalness of the LECs is realized but the 
systematicness is not. We should say that a low-momentum potential, like that 
represented by the dashed curve in Fig. |^ is simply a phenomenological model 
meant only to reproduce the scattering length and the effective range. Often, in 
the situation that a prediction made using the conventional NEFT begins to deviate 
from the data as the energy is increased, it is claimed that the energy is beyond the 
region to which NEFT is applicable. However, as we have seen, the reason that there 
is a small applicable energy region is that the conventional NEFT parametrization 
is not suitable to simulate the UT-generated model-space potential. 

Thirdly, the UT method does not preserve the values of the off-shell T-matrix 
elements. This leads to a (partial) breakdown of the EFT idea of the separation of 
scales. We now consider an example as an illustration. We start with the full space 
and consider a matrix element of an operator O corresponding to a two-nucleon state 
i.e., {^\0\^). First, we discuss the evaluation of this matrix element using the 
proper NEFT. In this method, we begin by decomposing the matrix element into 
low- and high-energy parts, as 

= + {^h\)0{\^l) + Wh)), (4-5) 

with 

= nW) , \^h) = A|iP') , (4-6) 

where rj and A are the projection operators defined in Eqs. (|2-7|) and (|2-8|) . respec- 
tively, and the subscripts L and H indicate the low- and high-energy components, 
respectively. Then, we proceed as 

{^\0\^) = {^l\0\^l) + the other terms (4-7) 

= {^l\0\^l) + {^Lpheavym 
~ {^l\0\^l) + {^L\Ocnt\^L) , 

where Oheavy is given by one insertion of O with all orders of rescattering terms due 
to the A^A'^-potential; all intermediate states involved are the high-momentum states 
contained in the projection operator A (see Fig.lJSlfor a diagrammatic representation 
of Oheavy)- Further, Ocnt is a series of contact operators and their couplings are 
generally dependent on the on-shell energies of the initial and final states. They 
are determined such that the best simulation of Oheavy is realized. The step from 

*-* The solid curve in Fig. |5|is obtained from the Viowk method. It is noted, however, that the 
Viowk method and the UT method yield very similar low-momentum potentials. 
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Fig. 23. A diagrammatic representation of Oheavy introduced in Eq. II4-7I I. An operator O acts 
on the full space, while Oheavy acts on the model space including only low-momentum states. 
The shaded blob includes the free two-nucleon propagator plus the ladder diagrams generated 
by the full-space A'^A''-potential, V. In all of the intermediate states involved in the diagrams, 
the nucleons run over the complementary space of the nucleon momentum states. 

the first to the second line in Eq. 1)4 • 7(1 indeed represents the model-space reduction 
according to the WRG equation. The step from the second to the third line is made 
by carrying out an expansion of Oheavy in terms of the contact operators. The third 
line is what we consider in NEFT. The low-energy physics, namely, that described 
by {^lIOI'Fl) is treated identically to the original theory. The separation of scales 
is clearly realized, which we expect for EFT. 

Next, we consider the case for the UT method. Here we have 

{^\0\^) = {^\UU^OUU^\^) (4-8) 
= {^'\{0 + V^O + OV + V^OV)\^') 

where U is the unitary operator introduced in Eq. (|'2-9() to reduce the model space, 
and we have used the definitions ^' = ^ rj^ and V = U — 1. The step 

from the first to the second line represents the model-space reduction using the UT 
method. The third line is what one evaluates in the conventional NEFT. As we 
see, the low-energy physics is not treated in the same way as in the original theory; 
the high-energy physics enters through {'I^'\0\'I^') , while the low-energy physics is 
mixed in (l^''|C'^„J'P''). In fact, the UT method causes mixing of scales rather than 
separation of scales. This is why, as stated above, a clear separation of scales cannot 
be accomplished with the UT method. Actually, the situation is even more serious. 
In the case that O is a Hamiltonian, as we have seen in our result, the step from the 
second to the third line in Eq. H4-8() (the expansion in terms of contact interactions) 
cannot be carried out satisfactorily. For the other operators, such as electromagnetic 
currents, the same difficulty may exist. 

The following arguments hold regardless of whether the conventional NEFT is 
related to the UT method. In the conventional NEFT, one starts from a model 
space. Thus, what one does in the conventional NEFT is not to follow all the steps 
taken in Eq. ()4-8() but to start from the third line. In this case, is generated by 
a potential whose role is only to reproduce on-shell quantities for a rather limited 
energy region, as discussed in the previous paragraph. This means that the off-shell 
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behavior of ^' is completely out of our control. What one can do in the conventional 
NEFT is to determine the couplings involved with O^nf ^ given kinematical point 
by fitting to data and hope that this uncontrollable behavior is completely canceled 
by Ocni- However, all the ad hoc steps taken in the conventional NEFT, as discussed 
above, are expected to shrink the kinematical region in which the theory is valid. 
The proper NEFT should have much broader applicability. 

Despite the weakness of the conventional NEFT described above, it has had 
much success. We now consider the reason it seems to have been successful. A likely 
reason is that most of the previous works have been concerned with EFT(7r), and a 
relatively large cutoff value was used.*^ As we have seen in Fig. [71**) the behavior of 
the UT-generated model-space potential around the cutoff becomes much smoother 
for a larger cutoff. Therefore, it is practically possible in this case to realize a more 
accurate simulation of the model-space potential in terms of one- (and two-) pion- 
exchange plus contact interactions. Furthermore, it is noted that Vm's obtained 
from the WRG and the Viowk methods are very similar for a larger cutoff, as we 
found in Fig. [7| Therefore, the conventional Veft{it) with a large cutoff can be 
interpreted as a parameterization of the WRG-based Vm, which also supports the 
practical use of the conventional procedure. On the other hand, it is not expected 
that VEFTifi-) constructed in the conventional manner is effective, because the pionless 
potential should be defined in a rather small model space. However, the usefulness 
of EFT should not depend on the value of the cutoff if appropriate d.o.f. for the 
cutoff are considered explicitly. Even though the conventional EFT(7r) with a large 
cutoff is practically useful, it is noted that this fact does not completely justify the 
conventional NEFT. 

Lastly, we make one more comment on the success of the conventional Veft- 
Although the scattering observables have been thoroughly studied with the use of 
Veft, electroweak processes have not,***) and the off-shell behavior of the wave 
functions generated by Veft have not yet been thoroughly examined. As we have 
discussed, even if some potentials can accurately predict the scattering phase shifts, 
they do not always control the off-shell behavior. We should consider the success of 
the conventional NEFT while keeping this point in mind. 

4.4. What have we learned? 

In this work, we have explicitly demonstrated that the RG idea is valid in ap- 
plication to the nuclear force. The validity of the RG idea is a condition for NEFT 
to be effective. The RG idea is that as we reduce the model space sufficiently, or 
equivalently, as we view a system in a sufficiently coarse-grained manner, we find 
the system to be controlled by a dynamics that are independent of the details of the 
small scale physics. Although we may need a confirmation, it is likely that the RG 



*' There is a well-known success of EFT(^) with KSW-counting,*' in which dimensional regu- 
larization plus the so-called power divergence subtraction is used. 

**' It is noted that the result for the UT method is essentially the same as that for the Viowk 
method. 

.**) Ai^jjQ^gjj there have been many NEFT-based studies of electroweak processes in few-nucleon 
systems, a Veft constructed in the conventional manner has not been used directly in those studies. 
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idea is valid in application to two-nucleon systems as a whole, i.e., in application 
to the nuclear current operator as well as the nuclear force. In the following, we 
consider what we can learn about a description of a two-nucleon system based on 
a traditional, phenomenological model approach, and about a description based on 
NEFT from a point of view rooted in the RG idea in a two-nucleon system. 

A model used in the traditional approach contains well-known large scale physics 
and phenomenological small scale physics. There are many such models which differ 
from each other with regard to small scale physics. Those models are able to describe 
a given low-energy reaction with similar precision. The RG approach tells us that 
these models are equivalent at low energies; all of these models are, after an appropri- 
ate model-space reduction, transformed into essentially a single low-energy theory, 
which does not have the dependence on the modeling of the small scale physics. 
This low-energy theory can be interpreted as NEFT. Therefore, the RG procedure 
guarantees that the models and NEFT give the same model-independent result in a 
low-cncrgy region in which the details of small scale physics involved in the models 
are unimportant. We do not have to concern ourselves with the model dependence 
of the predictions obtained with the models, as long as we are concerned with that 
energy region. This point of view based on the RG idea may be interpreted as a 
foundation of the traditional model approach. 

As we have seen, the RG idea makes the relation between NEFT and the model 
approach clear. Prom this viewpoint, we give a supplementary explanation for the 
conventional argument regarding the interplay between NEFT and the model ap- 
proach. Typically, models arc constructed in a model space that is much larger than 
those for NEFT. Therefore, these models are, in principle, applicable to energy re- 
gions outside the region to which NEFT is applicable. It has sometimes been claimed 
that it is a reasonable strategy to apply a model to a higher energy region if its re- 
liability had been tested in the low-energy region by comparing its predictions with 
NEFT-based ones. As we go into higher energy regions, however, the dependence on 
the modeling of small scale physics gradually enters into model predictions. Because 
the details of the small scale physics involved in the models has never been tested 
by the comparison, this type of extrapolation is expected to be rather limited. For 
a better extension of the model to higher energy regions, it is important to test the 
small scale physics part by using experimental data from that energy region. 

With the understanding of the relation between NEFT and the traditional model 
approach, we now discuss what is, and is not, newly gained by working with NEFT. 
One gain is that NEFT enables us to perform a perturbative calculation. As a 
consequence, we can improve the accuracy of a calculation systematically and we 
can estimate the theoretical uncertainty. It is noted, however, that there are cases 
in which the perturbation expansion does not converge sufficiently rapidly. This is 
because the contribution of the leading terms is suppressed by a certain symmetry. 
An example is the pp —>■ ppir^ reaction.^-*^) Another gain is that we can work with 
interactions free from the model dependence that results from treating the details 
of the small scale physics. This does not mean, however, that only NEFT gives a 
model-independent result. In the energy region in which NEFT is useful, NEFT 
and a reasonable phenomenological model are equivalent in the sense that they are 
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connected through the RG, as we have seen. ActuaUy, we should recognize that a 
new phenomenological prediction is rarely gained by working with NEFT. We can, 
in many cases, construct a reasonable model which is equivalent to NEFT at low 
energies; this is what nuclear physicists have done. Even so, there is one more gain, 
namely, that NEFT enables us to construct a nuclear force and nuclear electroweak 
current operator much more efficiently with much simpler parametrization than in 
the case that we construct a model. This efficiency is due to the fact that we know 
which interactions are to be considered at a given order of the NEFT perturbation. 
In this context, one may say that the chiral symmetry plays an important role 
in EFT(7r) because the chiral symmetry determines how the pion interacts with 
the other particles, and the use of interactions satisfying the chiral symmetry is 
important in setting up the counting rule.*) In constructing a model, on the other 
hand, one has to find an appropriate combination of interactions by trial and error, 
except for the well-known long-range mechanism. The simple parameterization of 
interactions in EFT results from the fact that we ignore details of small scale physics. 

§5. Summary and conclusion 

We have demonstrated the relation between Veft and Vph in this work. For 
this purpose, we studied how to reduce the model space. One method considered 
as a possibility for this reduction is the WRG equation, which is equivalent to the 
path integral method of integrating out heavier d.o.f. Wc found that the WRG 
equation generates a potential which can be accurately simulated by Veft] the 
basic assumptions of NEFT, namely, naturalness, systematicness and integrability 
are realized in Veft- Thus, we conclude that the WRG method is the appropriate 
model-space reduction in NEFT. Simultaneously, we conclude that Vph and Veft 
(and different Veft's acting on different model spaces) are connected through the 
WRG. It was shown that Veft is free from model dependence, i.e., dependence on 
the description used for the small scale physics. The use of simple contact interactions 
in NEFT can be naturally understood from the RG point of view. 

We also examined two other model-space reduction methods, the Viowk a-nd 
the UT methods. However, we found that they are not consistent with the path 
integral method and generate potentials that are not consistent with the basic NEFT 
assumptions. Wc conclude that they are not appropriate methods of model-space 
reduction in NEFT. 

*^ It may not be warranted to make this typical argument, from the resuh obtained in this work, 
about the role played by the chiral symmetry in NEFT. Our result implies that any framework (i.e., 
a Lagrangian and a perturbation scheme) is useful if it gives a good parameterization of a model- 
space A'^A'^-potential which is free from model dependence coming from the description of small scale 
physics; NEFT is one of them. In order to assess the role of the chiral symmetry, it would be 
helpful to examine an A'^A''-potential with A ~ 400 MeV in which the details of multi-pion-cxchango 
potentials are expected to play an important role. In this case, one has to explicitly include the 
multi-pion-exchange and therefore needs a rule to perturbatively include them; the use of the chiral 
Lagrangian may be important for this purpose. If the OPEP is the only important pion-induccd 
mechanism, as in our result, we do not have to rely on the chiral Lagrangian, because many NN- 
models include the OPEP without using the chiral Lagrangian. 
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After finding the relation between Veft and Vph-, we discussed the Veft proper- 
ties which have not yet been fuhy recognized. Because Veft is obtained by integrat- 
ing out heavy d.o.f. following the RG procedure, the Veft should have properties 
consistent with the model-space reduction. For example, Veft is on-shell energy 
dependent; we proposed a method for modifying the conventional counting rule. We 
emphasized the importance of being careful that Veft be defined in a model space 
whose size is sufficiently smaller than those used for Vph- We emphasized this point 
because there has been confusion that has resulted from overlooking this point. 

From the point of view developed in this work, we discussed the nature of the 
conventional Veft- We discussed the fact that the conventional Veft can be inter- 
preted as a parametrization of a potential whose model space has been reduced from 
some theory with the use of the UT method. Whether or not this interpretation is 
correct, the conventional Veft has problems from which the proper Veft is free. 
In particular, one serious problem is that the mixing of scales rather than the sepa- 
ration of scales is involved in the conventional Veft- This implies that a model is 
introduced into the theory to a certain extent. These problems appear clearly when 
the value of the cutoff is small. We conclude that the conventional Veft is not fully 
consistent with the basic idea of NEFT, but it is practically useful in the case of a 
large model space and has had been successful phenomenologically. 

Finally, we discussed what we can learn about NEFT and the traditional model 
approach from the RG idea. We conclude that the model approach is equivalent at 
low energies to NEFT in describing low-energy two-nucleon systems. Therefore, we 
do not have to be concerned with the model dependence of predictions made by a 
reasonable phenomenological model, as long as we are concerned with energy regions 
in which details of small scale physics is unimportant. We also conclude that we can 
gain the following by working with NEFT: a perturbative calculational procedure 
which enables us to systematically add higher-order corrections and to estimate 
theoretical uncertainty; a model-independent framework which, however, does not 
mean that only NEFT gives model-independent predictions; an efficient and simple 
construction of nuclear systems, for which chiral symmetry plays an important role 
in the case that the pion is dynamical. 
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Appendix A 

Derivation of the Wilsonian Renormalization Group Equation 



In this appendix, we derive the WRG equation [Eq. (|2-3|) ] by explicitly integrat- 
ing out the nucleon high-momentum states in a path integral. We restrict ourselves 
to the two-nucleon center of mass (CM) system. Therefore, we consider only a 
two-body force; that is, we do not consider the inclusion of intrinsic multi- (i.e., 
more than two) body forces nor their generation due to integrating out heavy d.o.f. 
For simplicity, we consider the case in which there is no mixing of different partial 
waves. Extension to the coupled-channel case is straightforward. We start with a 
path integral Z written only in terms of the nucleon field A^": 

Z = j VNVN^e'^'^^ , (A-1) 

with 

)(a)t TT/(a) n(") 



Cn = Co-Y, Dr ^^"^ , (A-2) 



a, A 

= (i9o + (A-3) 

where a specifies a partial wave of an A^A^-state and P^"^ denotes the projection 

operator onto a partial wave a. [Explicit expressions for P_J"^ are given in Eq. (A4) 
of Ref. 22).] The suffix A indicates the direction of the spatial and isospin-spatial 
polarization of the two-nucleon system. The quantity W^"^ is the coupling of the two- 
body A^A^-interaction for a partial wave a. W^°'^ contains derivatives and therefore 
depends on the magnitudes of the momenta carried by the free incoming nucleons (/c) 
and outgoing nucleons {k') at vertex; it can also depend on the on-shell momentum 
(p) of the two-nucleon system and the momentum cutoff (A). We can thus write 
= W^^'^k' ,k;p, A). The arguments of are suppressed until Eq. HA-12|) 

for simplicity. The quantity V^"'^ defined by V^"^ = W^"^ /2 is interpreted as an 
A^A^-potential. 

We separate the nucleon field into high-frequency modes Nh and low-frequency 
modes A'^^. We then expand the high-frequency modes of the interaction Lagrangian 
to obtain 

Z = J VNlVNI e' ^ '^^ ^ J VNhVNIj e' ^ '^^ ^^"^ 1 1 (A-5) 
{-if 



+ 



2! 

a/3,Ap 
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where includes only the nucleon fields with low- frequency modes. Integration 
of the first term in Eq. HA-5|) over Nh leads to an unimportant overall factor of 
Z. We explicitly perform the integration of the second term in the following. This 
integration leads to the one-loop correction to the two-body interaction. The ellipsis 
includes all possible diagrams of products of connected multi-loop ladder pieces. The 
self-energy pieces are not considered, because we should maintain consistency with 
the NEFT calculational method. In the NEFT calculation, one solves the Lippmann- 
Schwinger equation, which means that only the ladders of the A^A^-potential are 
resummed. 

Integrating the second term above yields 

j VNhVN\j e^/'^^ ^''"^ dP^{xi)D^j^'^^ {X2) j j VNhVN\j e'^'^'^ ^"""^ (A-6) 



H (27r)4 (27r)4 IgO - ql/2M + ieq^- ql/2M + ie 



2Tr ( P,^" 



A 



where indicates that qi [i = 1,2), the nucleon momentum in the intermediate 
state, runs over the region for the high-frequency modes. In Eq. (|A-6() . we have used 

the relation = -P'^^ Now, Eq. (U^ becomes, up to an unimportant overall 

factor, 

2Tr fp^^^P^^^''"! e''^ii+i2)<x2-xi) 



{q'l- qI/2M + ie) {q^ - ql/2M + ie) "^^^ " 



I VNlVNI e'l'''^ ^N^li + i(27r)45(4)(P - P') Z)J"/(0)Ty(") 



d^, 2Tr (Pf)pr^ 



H (2vr)4 (gO _ q2j2M + te) (pO - qO - (P - q^)^ /2M + ie) 
j VNlVNI e*/'^'- - ^(27r)V4)(P - P') 

2^ ^LX iU)>^ (27r)3 pO - qyu + ie^ + 



|^(/3)z)|/^)(0) + ... 



a/3,Xp 



where we have used the fact that we are working in the CM system. The energy and 
momentum of the initial (final) two- nucleon system is denoted P = {P^,P) (P')- 
We consider the case in which an infinitesimally small momentum shell is integrated 
out; that is, qi in Eq. ()A-7|) runs over a region A — 5A < \qi\ < A. In this case, use 
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of the relation / dQq^2Ti (^p|°^Pp^^^^) = AnSap 5xp leads to 

Z = j VNlVnIc'-^'^^ '^N^\^-l^{2TTf5^^\P-P') (A- 



with 



47r2 PO-^VM ^ ' 



By resumming the terms in the curly brackets in Eq. ()A-8p as an exponential, we 
obtain the path integral 

Z = j VNlVNI e'^'^"'' ^^^'-ff , (A- 10) 



with 

^N,eff 

= - i D^fix) + + O {{6Af)) 



-C^i^-iJ2DtIH^)(sw(-^+0{{5A)'))Di%) (A-ll) 



(a) 



■l X y^i 



Therefore, for an infinitesimally small reduction of the momentum cutoff, we have 
the following renormalization group equation for W^"^: 

Here, the minus sign on the l.h.s. indicates that the r.h.s. is the shift of VF^") due 
to an infinitesimal decrease of A. We have also used = p'^/M. Thus we obtain 
the WRG equation in terms of introduced above: 

Appendix B 

NN -Potential 



We now present expressions for Veft{tt) used in this work in the partial- wave 
basis. These are obtained from Eq. H2-14() by changing the basis. The expres- 
sions presented here are substituted directly into the Lippmann-Schwinger equation, 
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Eq. (E31). The result for the channel, F^'^o) = (^^ol^l^^o), is 

+ Ci'^°) + C^'''°\k^ + k'^) + C7f {{k^ + k'^f + {2kk'f/^] , 

with 

Lj{^)^ j\t^^ , (B-2) 

and 7 = (A;2 + A;'^ + ml)/2kk'. Here, Pj(t) is a Legendre function of order J. The 
other symbols are the same as those used in the text. We use the average value for 
the pion mass, m^r = (m'7r+ + "i^- + mjro)/3, and we adopt the phenomenological 
value g\/Af'^ = 0.075 • 47r/m^. The LECs of the contact interactions depend on 
A and p. The values of the LECs for some sets of A and p are given in Table Q 
The result for ^^\k' ,k;p, A) = {k',^Di\V\k,^Si) is obtained by simply replacing 
Cf in Eq. dEU) with cf {i = 0,2,4). The result for yC^-^i) = {^Di\V\^Si) 
is 

V('''-'^\k',k;p,A) = ||^/2 |-Li(7) + ^^0(7) + ^^2(7)} (B-3) 

+ D^^'h'^ + D[''h'^ l^h^ + k'^^ , 

and interchanging k and k' on the r.h.s. gives the expression for {k' ,^ Si\V\k,^ Di) . 
Finally, the result for F^'^i) = {^Di\V\^Di) is 
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